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We derive several results concerning non-perturbative renormalization in the 
spherical field formalism. Using a small set of local counterterms, we are able 
t— i , to remove all ultraviolet divergences in a manner such that the renormalized 

theory is finite and translationally invariant. As an explicit example we 
consider massless cp 4 theory in four dimensions. [PACS numbers: ll.10.-z, 
ll.10.Gh, 11.15.Tk] 
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1 Introduction 

Spherical field theory is a non-perturbative method which uses the spheri- 
cal partial wave expansion to reduce a general rf-dimensional Euclidean field 
j> ! theory into a set of coupled radial systems ([0], [fj). High spin partial waves 

correspond with large tangential momenta and can be neglected if the theory 
is properly renormalized. The remaining system can then be converted into 
differential equations and solved using standard numerical methods. <^ 4 the- 
ory in two dimensions was considered in . In that case there was only one 
divergent diagram, and it could be completely removed by normal ordering. 
In general any super-renormalizable theory can be renormalized by removing 
the divergent parts of divergent diagrams. Using a high-spin cutoff J max and 
discarding partial waves with spin greater than J max , we simply compute the 
relevant counterterms using spherical Feynman rules. 
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The J m ax cutoff scheme however is not translationally invariant. It pre- 
serves rotational invariance but regulates ultraviolet processes differently de- 
pending on radial distance. In the two-dimensional (p A example it was found 
that the mass counterterm had the form 

Cct. oc 2 (t) [K {flt)l (flt) + 2E n= i, Jmax ^n(^)4(^)] , (1) 

where I n , K n are n th -order modified Bessel functions of the first and second 
kinds, [i is the bare mass, and t is the magnitude of t. As J max — > oo, we find 

C c ,. oc m [log(^) + 0(J m D] . (2) 

Our regularization scheme varies with t, and we see that the counterterm also 
depends on t. The physically relevant issue, however, is whether or not the 
renormalized theory is independent of t. In this case the answer is yes. Any 
t dependence in renormalized amplitudes is suppressed by powers of ./"ax, 
and translational invariance becomes exact as J max — > oo. 

We now consider general renormalizable theories, in particular those which 
are not super-renormalizable. In this case the number of divergent diagrams 
is infinite. Since we are primarily interested in non-perturbative phenom- 
ena, a diagram by diagram subtraction method is not useful. In the same 
manner strictly perturbative methods such as dimensional regularization are 
not relevant either. Our interest is in non-perturbative renormalization, 
where coefficients of renormalization counterterms are determined by non- 
perturbative computations.^] In this paper we analyze the general theory 
of non-perturbative renormalization in the spherical field formalism. In the 
course of our analysis we answer the following three questions: (i) Can ultra- 
violet divergences be cancelled by a finite number of local counterterms? (ii) 
Can the renormalized theory be made translationally invariant? (iii) What 
is the general form of the counterterms? 

The organization of the paper is as follows. We begin with a discussion of 
differential renormalization, a regularization-independent method which will 
allow us to construct local counterterms. Next we describe a regularization 
procedure which is convenient for spherical field theory. In the large radius 
limit t — > oo our regularization procedure (which we call angle smearing) is 

We should mention that Pauli-Villars regularization is compatible with non- 
perturbative renormalization. However this introduces additional unphysical degrees of 
freedom and tends to be computationally inefficient. 
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anisotropic but locally invariant under translations. For general t we expand 
in powers of t^ 1 to generate the general form of the counterterms. We con- 
clude with two examples of one-loop divergent diagrams. We show by direct 
calculation that the predicted counterterms render these processes finite and 
translationally invariant. 



2 Differential renormalization 

Differential renormalization is the coordinate space version of the BPHZ 
method.f] It is framed entirely in coordinate space, and renormalized am- 
plitudes can be defined as distributions without reference to any specific 
regularization procedure. Differential renormalization was introduced in 0L 
and a systematic analysis of differential renormalization to all orders in per- 
turbation theory using Bogoliubov's recursion formula was first described 
in 0. The usual implementation of differential renormalization is carried 
out using singular Poisson equations and their explicit solutions. In our 
discussion, however, we find it more convenient to operate directly on the 
distributions.^ We describe the details of our approach in the following. We 
should stress that the two approaches are equivalent, differing only at the 
level of formalism. 

We assume that we are working with a renormalizable theory. For indices 
ii , ■ ■ ■ ij let us define 

f 1 --^- =t h t i2 (3) 
V n ,... b . = V n V 42 ---V, r (4) 

Let fit) be a smooth test function, and let I(t — t 1 ; jj 2 ) be a smooth function 
with support on a region of scale yu -1 . We define Si, [/] (t) as I(t — i 1 ;^ 2 ) 
multiplied by the j th order term in the Taylor series of f(i) about the point 
Inserting delta functions, we have 

SU\t) = I(t- ? ; /i 2 ) E V^^V^Ji?)} (5) 



z . 



°Paraphrase of private communication with Jose Latorre. 

6 Our approach is similar to the natural renormalization scheme described in [^). In 
contrast with p], however, we do not a priori specify the finite parts of amplitudes. 
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For the purposes of this discussion we will require 

I(t- ?; n 2 ) = 1 + N {t- ?) as ? -> t, (6) 

where iV is some positive integer greater than the superficial degree of di- 
vergence of any subdiagramQ in the theory we are considering. For any 
renormalizable theory N > 2 will suffice. In our formalism, lit — r ; fx 2 ) de- 
termines how finite parts of renormalized amplitudes are assigned, and fi is 
the renormalization mass scale. 

We now consider a particular diagram, G, with n vertices. We define 
K(ti, ■ ■ ■ t n ) to be the kernel of the amputated diagram, i.e., the value of the 
diagram with vertices fixed at points ti, ■ ■ ■ t n . The amplitude is obtained by 
integrating K(t\, ■ ■ -t n ) with respect to all internal vertices. We will regard 
K as a distribution acting on n smooth test functions f±, ■ ■ • f n . (For external 
vertices containing more than one external line and/or derivatives, f e xt{t e xt) 
should be regarded as a product of test functions, with possible derivatives, 
at t ex t.) 

K : fi, ■ ■ ■ f n -> Jd%. • • • d% K(h, ■ ■ ■ f n )/i(ti) ■ ■ • f n (t n ). (7) 

Let us assume that our diagram is primitively divergent with superficial de- 
gree of divergence j. We now define another distribution TqK, which extracts 
the divergent part of K. We start with the case when G has more than one 
vertex. Let us define TqK : fi, ■ ■ ■ f n — > 



E Jd%---d 4 t n K(t 1 ,---t n )Sp fi(h)'.-Sp f n (t n ), (8) 

Lave 



31 H Vjn<3 



where t ave = -{t\ + • • ■ + t n ). We note that the subtracted distribution K — 
TqK is finite and well-defined for all fi, ■ ■ ■ f n . Let us define 

= Jd% ■ ■ ■ d% 5 4 (fc±S, _ T)K{t 1: ■ ■ ■ t r: 



n 

k=l,--n 



i(t k -t-y)(t k -t) ^ 



7 In our discussion a subdiagram is a subset of vertices together with all lines contained 
in those vertices. 
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We can then rewrite T G K : fi,---f n 



E E 

Jl + ---+jn<j 11,1,42,1'" 
rf 'l,n'"' l 'j n ,n 



JdH^ 



»l,lj*2,l - "» 

K 



n 



k=l,---n ' H,fc>—* 



(t) J <i 4 Fi ■ • • (i 4 z n 
k S 4 (t- z k )J • • • / n (4) 



(10) 



The delta functions make this kernel completely local. We can read off the 
corresponding counterterm interaction by functional differentiation with re- 
spect to each of the component functions of f ex t(t ex t) for the external vertices 
and setting fmtitint) = 1 for the internal vertices. We now turn to the case 
when G has only one vertex. For this case we set TqK = K, which is equiva- 
lent to normal ordering the interactions in our theory. In this case K is itself 
local and therefore TqK and our counterterm interaction are again local. 

We now extend the definition of Tq in (ID) to include the case of subdia- 
grams. Let G be a general 1PI diagram, and let G' be a renormalization part| 
of G with superficial degree of divergence f . For notational convenience we 
will label the vertices of G so that the first n' vertices lie in G' . If G' has 
only one vertex then again we normal order the interaction. Otherwise we 
define T G ,K : f 1} ■ ■ ■ f n -> 



E fd%---d%K(t u ---t r 



3'i+-+3',<3' 



tav 



fn'+l{tn'+l 



■ •si" 



fn' (t n > 

e 

fn (^n) 



where t ave = ^(tx + • ■ ■ + t n /).f\ This definition can be used recursively to 
define products of Tq>Tq> 2 for disjoint subdiagrams G' x D G' 2 = or nested 
subdiagrams G[ D G' 2 . For the case of nested subdiagrams we always order 
the product so that larger diagrams are on the left. 

It is straightforward to show that the T operation acts as the identity on 
local interactions and thus treats overlapping divergences in the same manner 
as BPHZ. Following the standard BPHZ procedure (|§ — f|), we can write 
Bogoliubov's R operation using Zimmerman's forest formula, 



R 



e n (-V; 

F 7SF 



;i2) 



8 A renormalization part is a 1PI subdiagram with degree of divergence > 0. 
9 After applying Tg', we regard G as being contracted to single vertex at t a , 



5 



where F ranges over all forests^ of G, and 7 ranges over all renormalization 
parts of F. In the product we have again ordered nested subdiagrams so that 
larger diagrams are on the left. Let j be the superficial degree of divergence 
of G. The renormalized kernel, RK, is given by 

RK = RK forj<0 , , 

RK ={l- T G )RK for j > 0. ^ ' 

Our final result is that all required counterterms are local, and the form of 
the counterterms is 

£c, = £^,v^) Fa®M<K%> V 4 0(t5), (14) 

where the sum is over operators of renormalizable type. For the case of gauge 
theories, our renormalization procedure is supplemented by the additional re- 
quirement that the renormalized amplitudes satisfy Ward identities.^ If our 
regularization procedure breaks gauge invariance these identities are not au- 
tomatic and the required local counterterms will in general be any operators 
of renormalizable type (not merely gauge-invariant operators). This is, how- 
ever, a separate discussion, and the details of implementing Ward identity 
constraints will be discussed in future work. 



3 Regularization by angle smearing 

In this section we determine the functional form of the coefficients Fa(c) i n 
flnp. To make the discussion concrete, we will illustrate using the example 
of massless 4 theory in four dimensions 

C=\4>V 2 4>-^ + C c . t .. (15) 
From (jrjl) £ c t is given by 

F 02 (t)0 2 (t) + E«4v#(*1 V '^Wl + ( 16 ) 

Let G(t, v ) be the free two-point correlator. We will use a regularization 
scheme which preserves rotational invariance and is convenient for spherical 

10 A forest is any set of non-overlapping renormalization parts. 

11 See (To) for a discussion of gauge theories using the method of differential renor- 
malization. 
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field theory, but one which breaks translational invariance. We regulate 
the short distance behavior of G by smearing the endpoints over a radius 
t spherical shell within a conical region R,M 2 (i), where R,M 2 (t) is the set of 
vectors u such that the angle between t and u is between — and (see 
Figure 1). The result is a regulated correlator 

G M 2(t,i > ) = 1 ^p-i / n ^ d 3 ud 3 u'G(tu,t'u'). (17) 

u'eR M 2(?) 

We recall that our renormalized theory is determined by the translationally 
invariant function /(t — t ave ; /i 2 ) described in the previous section. Even 
though our regularization scheme breaks translational invariance, the renor- 
malized theory nevertheless remains invariant. 

As the radius t increases the curvature of the angle-smearing region be- 
comes negligible. In the limit t — > oo the region becomes a flat three- 
dimensional ball with radius lying in the plane perpendicular to the radial 
vector. In this limit our regularization is invariant under local transforma- 
tions and the counterterms converge to constants independent of t, 

I^F^i) = 4%(^) (18) 
lim^® = A^$>(£) (19) 

t— >oo ^ 

ImF^® = #>(£). (20) 

We have chosen our coefficients to be dimensionless. Although our reg- 
ularization scheme is invariant under rotations about the origin, the radial 
vector has a special orientation which is normal to our three-dimensional ball. 
Our regularization scheme is therefore not isotropic. The result (as should 
be familiar from studies of anisotropic lattices) is that the coefficient of the 
kinetic term has two independent components 

ij,(o) _ tt,(o) & (o) , ^ (21) 

Starting with the t —> oo result at lowest order, we now expand our 
coefficient functions in powers of 

pij (f\ — r vM ( M 2 A _|_ J^J.C 1 ) fl» 2 1, 1 r V,( 2 ) { \ _i foo\ 

F, 2 (t) = M^(^) + f + + • • • (23) 

= + M^(^) + W c ^(^) + • • • • ( 24 ) 
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For the moment let us assume t > A 1 for 

A = mgM 1 " 2 , (25) 

for some fixed mass m and constant z such that < z < \. In this region 
our dimensionless expansion parameter is bounded by ^| and therefore 
diminishes uniformly as M — >• oo. 

In general the -jrp dependence in the functions will contain analytic 
terms as [i 2 —> as well as logarithmically divergent terms. There are, 
however, no inverse powers of These would indicate severe infrared 

divergences not present in the processes we are considering, as can be deduced 
from the long distance behavior of the integral in With this we can 

neglect terms which vanish as M — > oo, 

Fp®=M*c<$(&) + (26) 

F*® = (28) 

Since our regularization scheme is invariant under M — > —M, we have also 
omitted the term proportional to c? 2 which is odd in M. 

We now consider what occurs in the small region near the origin, t < A -1 . 
For the theory we are considering (and in fact for any renormalizable theory) 
the highest ultraviolet divergence possible is quadratic]^ In the limit M — > 
oo we deduce that each Fa scales no greater than 0(M 2 ). On the other hand 
the volume of the region t < A -1 diminishes as 0(M 4z ~ 4 ). Thus the total 
contribution from the region t < A -1 scales as 0(M 4z ~ 2 ) and can be entirely 
neglected. 

To summarize our results, the counterterm Lagrange density has the form 



4°U(v0(t)) 2 + 4fi(i ■ v0(t)) 2 + {M 2 cf + M 2 2 V(t) + #y (*) 



(29) 



12 If our theory contained bare masses rrn, similar arguments would apply for the infrared 

2 

limit /i 2 ,?7i? — > 0, with — *h fixed. 

13 There may be additional logarithmic factors but this does not matter for our purposes 
here. 
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4 Spherical fields 



We now examine the results of the previous section in the context of spherical 
field theory. We start with the spherical partial wave expansion, 



<t> = E«=0,l,-^=0,-/E m =-n,...A™,m(t)^,n,m(^,^,^), 

where Yi jm>n are four- dimensional spherical harmonics satisfying 

/ d 3 Q Yy n , m , (9, ijj, if)Yi tTljm (9, ij),ip) = 8i, t i5 n , jn 5 m , tm , 



(30) 



(31) 



Y*{9^^) = {-l) m Y l ^ m {9^^). 



(32) 



The explicit form of l^ m ,n can be found in The integral of the free 

massless Lagrange density in terms of spherical fields is 



jd A tC = J °°dt 




H,n,m 



(33) 



It can be shown that the process of angle smearing the field <j>(t) is equivalent 
to multiplying (f>i )n>m (t) by an extra factor sf^t) where 



-M 



(*) 



2Mt[(Z+2)sin( I i I )-Zsin(af)] 
*(J+l)(Z+2)|2-Mtsm( 



(34) 



For large /, sf*(t) diminishes as l~ 2 , and so the correlator receives an extra 
suppression of /" 4 . We will later use this result to estimate the contribu- 
tion of high spin partial waves. The regularization of our correlator can be 
implemented in our Lagrange density by dividing factors of s z M (t), 



at 2 at 2 l \ b ^ 



'Z.n.m 



(35) 



We now include the interaction and counterterms. We first define 



h,ni,mi;l2,n2,m2 
I3 ,ri3 ,m3 ;/4 ,n4 ,wi4 

/3 



(36) 



Jd ttY hjnumi (6, ip, (p)Y htn2tm2 (9, ip, (p)Y h>n3>m3 (9, i/j, <p)Y h , n4}7n4 (9, i/j, <p). 



14 [[ll] deserves credit as the first discussion of radial (or covariant Euclidean) quantiza- 
tion, an important part of the spherical field formalism. 
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We can write the full functional integral as 



P0exp [jd'tC] OC J (lll,n,m^l,n,m) <*p [/ °°<ft + L 2 + L 3 )] , (37) 

where 



m-|Z(Z + 2) [(^(t))" 1 ^] 



l,m,n 



(38) 



E 



-i)' 



Ln,—m 



—c 



(0) 



— c 



tt,(0) 
V0V</> 



at 2 

1.2 



,2 ; 



Ln,m 



(39) 



3 I41 04 J / ,, [i3,m3,n3;i4,m,4,? 



,mi ,ni ^2 ,ra 2 ,n2 ( t ) lz ,m 3 ,n 3 ,m 4 ,n 4 ■ 



We have used primes in preparation for redefining the fields, 



sl,n,m- 



(40) 



(41) 



The Jacobian of this transformation is a constant (although infinite) and 
can be absorbed into the normalization of the functional integral. Now the 
Lagrangian L\ has the usual free-field form in terms of 0z, n , m while L2 and 
L 3 are now functions of sf 1 (i)</>z, n ,m- 

With M serving as our ultraviolet regulator, the contribution of high-spin 
partial waves decouples for sufficiently large spin I. We can estimate the order 
of magnitude of this contribution in the following manner. We first identify 
t L l (where t is the characteristic radius we are considering) as an estimate 
of the magnitude of the tangential momentum, pt- For p? 3> M ^> t^ 1 our 
correlator scales as By dimensional analysis, a diagram with N L loops 
and Nj internal lines will receive a contribution from partial waves with spin 
> I of order 



M 4 



(Pt) 



4Nr 



M 4 



(r'O 



1 A 4iVL 



(42) 
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5 One-loop examples 



We will devote the remainder of our discussion to computing one-loop spher- 
ical Feynman diagrams as a check of our results. Our calculations are done 
both numerically and analytically. The diagrams we will consider are shown 
in Figures 2 and 3. We start with the two-point function in Figure 2. The 
amplitude can be written as t 3 B{t) where 

B(t) oc E^^(^W) 2 - (43) 

Constants of proportionality are not important here and so we will define 
B{t) to be equal to the right side of ([43]). Our results tell us that if we 
choose our mass counterterms appropriately, the combination 

B(t) + M 2 ^ + (44) 

should be independent of t, or more succinctly, 

B(t) + (45) 

is independent of t. Let us first check this analytically. In the absence of a 
high-spin cutoff, we can explicitly calculate the sum in fl43|) : 

B(t) = £ + b{t) (46) 

where 

°W - ( 2 -Mtsin(^))2- W 

In the limit M — > oo, 

B{t)^^ + \M\ (48) 

We conclude that c^ 2 = —1 and B(t) + p-c^ 2 is in fact translationally invari- 
ant. 

In Figure 4 we have plotted B(i) — 4, computed numerically for various 
values of the high-spin cutoff J max - We have also plotted the limiting values 
b{t) and |M 2 . In our plot t is measured in units of m" 1 and B(t) — 4 is 
in units of m 2 , where m is an arbitrary mass scale such that M = 3m. As 
expected, the errors are of size There is clearly a deviation from |M 2 

''max ^ 

for t < M" 1 but the integral of the deviation is negligible as M — » oo. 
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We now turn to the four-point function in Figure 3. The amplitude can 
be written as t\t\C{ti,t 2 ) where 

C(ti, h) oc Ez,n, m (S " (tl (Sf (t2))2 IjhOih - h) + Ae(h - h)] 2 . (49) 



Again constants of proportionally are not important and so we will define 
C(ti, t 2 ) to be equal to the right side of (|4~9|). We can write C(t 1 , t 2 ) in terms 
of the regulated correlator G M 2(ti, t 2 ),f^ 

C(t u t 2 ) oc J d%d% [G M 2(t u t 2 )] 2 oc J d% [G M *{txM)\ 2 ■ (50) 

Since the coupling constant counterterm 

c$5%-t 2 ) (51) 

is translationally invariant, the amplitude by itself should be translationally 
invariant. Let us define 

J d%e~'^-^ [Gm^M)] 2 = f(p), (52) 

so that 

J d%e^ [G M i{t u t 2 )\ 2 = (53) 
Integrating over t\, we find 

J dhtZJiMCfah) oc ^Jxip^fif). (54) 

Let us define 

C{t) = $dt 2 tlJ x {pt 2 )C(t,t 2 )- (55) 
We now check that in fact 

C(t) oc ^Jxiph). (56) 



15 We recall that the regulated correlator goes with 4>'i nm rather than <j>i t n,m- But this 
is not important here since cf>' = </>o,o,o- 
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In the absence of a high-spin cutoff, we find that C(t) is given byQ 



C(t) = i J^) [ilog^ + cj +■•• , 
where the ellipsis represents terms which vanish as M — > oo and 

'1/2 



324 



17 / (sin k — k cos k) 
4fc 13 



324fc / ~ 



1/2 



7T (sin k—k cos fc) 



(57) 



(58) 



In Figure 5 we plot C(t) for different values of the high-spin cutoff J max as 
well as the large-M limit value 



d(t) = iJi(pti) Jlog^ + c 



(59) 



In our plot £ is measured in units of p 1 and M = 3p. From (|42|) the 
expected error is of size 



A/ 8 t 8 



We see that the data is consistent with the 



results expected. Again the deviation for t < M 1 integrates to a negligible 
contribution as M —>■ oo. 



6 Summary 



We have examined several important features of non-perturbative renormal- 
ization in the spherical field formalism and answered the three questions 
posed in the introduction. Ultraviolet divergences can be cancelled by a fi- 
nite number of local counterterms in a manner such that the renormalized 
theory is translationally invariant. Using angle-smearing regularization we 
find that the counterterms for 4 theory in four dimensions can be parame- 



terized by five unknown constants as shown in (|29|). Aside from our remarks 
about Ward identity constraints in gauge theories, the extension to other 
field theories is straightforward. We hope that these results will be useful for 
future studies of general renormalizable theories by spherical field techniques. 
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16 This calculation is somewhat lengthy. Details can be obtained upon request from the 
authors. 
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Figure 1. Sketch of the angle-smearing region (three-dimensional rendering). 
Figure 2. One-loop two-point correlator for </>o,o,o- 
Figure 3. One-loop four-point correlator for 0o,o,o- 
Figure 4. Plot of B{t) - 
Figure 5. Plot of C(t). 



14 



Figure 1 




Figure 2 



(Pl,n,m 




t 



Figure 3 




15 



Figure 4 




16 



